arXiv:1507.04780v2 [math.OC] 22Jul2015 


Distributed Average Tracking for Double-integrator Multi-agent 
Systems with Reduced Requirement on Velocity Measurements 


Sheida Ghapani ^ , Wei Ren ^ , Fei Chen , Yongduan Song 

’^Department of Electrical and Computer Engineering, University of California, Riverside, CA 92521, USA 
^Department of Automation, Xiamen University, Xiamen, China, 361005 
’’School of Automation, Chongqing University, Chongqing, China, 400044 


Abstract 

This paper addresses distributed average tracking for a group of physical double-integrator agents under an undirected graph with reduced 
requirement on velocity measurements. The idea is that multiple agents track the average of multiple time-varying input signals, each of 
which is available to only one agent, under local interaction with neighbors. We consider two cases. First, a distributed discontinuous 
algorithm and filter are proposed , where each agent needs the relative positions between itself and its neighbors and its neighbors’ 
filter outputs obtained through communication but the requirement for either absolute or relative velocity measurements is removed. The 
agents’ positions and velocities must be initialized correctly, but the algorithm can deal with a wide class of input signals with bounded 
acceleration deviations. Second, a distributed discontinuous algorithm and filter are proposed to remove the requirement for communication 
and accurate initialization. Flere each agent needs to measure the relative position between itself and its neighbors and its own velocity 
but the requirement for relative velocity measurements between itself and its neighbors is removed. The algorithm can deal with the case 
where the input signals and their velocities and accelerations are all bounded. Numerical simulations are also presented to illustrate the 
theoretical results. 
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1 INTRODUCTION 

This paper studies the following distributed average track¬ 
ing problem: given a group of agents and one time-varying 
input signal per each agent, design a control law for the 
agents based on local information such that all the agents 
will finally track the average of these input signals. The prob¬ 
lem has found applications in distributed sensor fusion [1], 
feature-based map merging [2], distributed Kalman filter¬ 
ing [3], where distributed computation of multiple time- 
varying signals are involved. Compared with the consensus 
problem, distributed average tracking poses more theoreti¬ 
cal challenges, since the tracking objective is time-varying 
and is not available to any agent. 

In the literature, linear distributed algorithms have been 
employed for special kinds of time-varying input signals. 
Ref. [4] uses frequency domain analysis to study consensus 
on the average of multiple input signals with steady-state 
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values. In [5], a proportional algorithm and a proportional- 
integral algorithm are proposed to achieve distributed aver¬ 
age tracking with a bounded tracking error, where accurate 
estimator initialization is relaxed in the proportional-integral 
algorithm. In [6], the internal model principle is employed 
to extend the proportional-integral algorithm to a special 
group of time-varying input signals with a common denom¬ 
inator in their Laplace transforms, where the denominator 
also needs to be used in the estimator design. Ref. [7] ad¬ 
dresses discrete-time distributed average tracking of time- 
varying input signals whose nth order difference is bounded 
with a bounded error. In [8], the authors propose Ist-order- 
input and 2nd-order-input consensus algorithms to allow the 
agents to track the average of their dynamic input signals 
with a pre-specified rate, where the interaction is described 
by a strongly connected and weight-balanced directed graph. 
In [9], a distributed continuous average tracking problem is 
addressed with some steady-state errors in , where the pri¬ 
vacy of each agent’s input signal is preserved. 

However, linear algorithms cannot ensure distributed aver¬ 
age tracking for general input signals. Therefore, some re¬ 
searchers employ nonlinear tracking algorithms. In [10], a 



class of nonlinear algorithms is proposed for input signals 
with bounded deviations, where the tracking error is proved 
to be bounded. A nonsmooth algorithm is proposed in [11], 
which is able to track arbitrary time-varying input signals 
with bounded derivatives. All the above references primar¬ 
ily study the distributed average tracking problem from a 
distributed estimation perspective, where the agents imple¬ 
ment local estimators through communication with neigh¬ 
bors freely without the need for obeying certain physical 
agent dynamics. However, there are applications where the 
distributed average tracking problem is relevant for design¬ 
ing distributed control laws for physical agents. One applica¬ 
tion is the region-following formation control [12], where a 
swarm of robots is required to move inside a dynamic region 
while keeping a desired formation. Then the dynamics of 
the physical agents must be taken into account in the control 
law design and the dynamics themselves introduce further 
challenges to the distributed average tracking problem. For 
example, the control law designed for physical agents with 
single-integrator dynamics can no longer be used directly 
for physical agents subject to more complicated dynamic 
equations. Distributed average tracking for physical agents 
with double-integrator dynamics is studied in [13], where 
the input signals are allowed to have bounded accelerations. 
Distributed average tracking for physical agents with gen¬ 
eral linear dynamics is addressed in [14] and [15]. Ref. [15] 
proposes a discontinuous algorithm, while a continuous al¬ 
gorithm is employed in [14] with, respectively, static and 
adaptive coupling strengths. Also in [16] a proportional- 
integral control scheme is extended to achieve distributed 
average tracking for physical Euler-Lagrange systems for 
two different kinds of input signals with steady states and 
with bounded derivatives. 

It is noted that in [13] both relative position and relative 
velocity measurements are required in the control laws de¬ 
signed for double-integrator agents. While double-integrator 
dynamics can be viewed as a special case of general lin¬ 
ear dynamics, the distributed average tracking algorithms 
in [15,16], when applied to double-integrator systems, still 
need relative position and velocity measurements. In prac¬ 
tice, velocity measurements are usually less accurate and 
more expensive than position measurements. In addition, 
relative velocity measurements are often more challenging 
and expensive than absolute velocity measurements. We are 
hence motivated to solve the distributed average tracking 
problem for physical double-integrator agents with reduced 
requirement on velocity measurements. This paper expands 
on our preliminary work reported in [?]. In the context of dis¬ 
tributed average tracking, reducing velocity measurements 
poses significant theoretical challenges. The reason is that 
unlike the consensus or single-leader coordinated tracking 
problems, there are significant additional inherent challenges 
in distributed average tracking as none of the agents has the 
tracking objective available. 

The contribution of this paper is summarized as follows. 
Two distributed algorithms (controller design combined with 
filter design) are introduced to achieve distributed average 


tracking with reduced requirement on velocity measure¬ 
ments. Each algorithm has its own relative benefits and is 
feasible for different application scenarios. In the first algo¬ 
rithm design, there is no need for either absolute or relative 
velocity measurements. Each agent’s algorithm only em¬ 
ploys its local relative positions with respect to neighbors, 
its neighbors’ filter outputs accessed through communica¬ 
tion, and the acceleration of its own input signal. The algo¬ 
rithm allows the agents to track the average of a large class 
of time-varying input signals with bounded acceleration de¬ 
viations, provided that the agents are correctly initialized. 
In the second algorithm design, there is no need for correct 
initialization and relative velocity measurements. Also inter¬ 
agent communication is not necessary and the algorithm can 
be implemented using only local sensing, which is desirable 
for certain applications (e.g., deep-space spacecraft forma¬ 
tion flying) where communication might not be desirable or 
available. Each agent’s algorithm only employs its local rel¬ 
ative positions with respect to neighbors, its own velocity, 
and its own input signal. Distributed average tracking can be 
achieved provided that the input signals and their velocities 
and accelerations are all bounded. 

Notations: Throughout the paper, M denotes the set of all 
real numbers and IR+ the set of all positive real numbers. 
Let 1„ and 0„ denote the n x 1 column vector of all 
ones and all zeros respectively. Let Aniax(') and Ainin(') 
denote, respectively, the maximum and minimum eigenval¬ 
ues of a square real matrix with real eigenvalues. We use 
0 to denote the Kronecker product, and sgn(-) to denote 
the signum function defined componentwise. Eor a vector 
function x{t) : M i—>■ M™, define ||a;(f)||p as the p-norm, 
x{t) e L 2 if x{t)'^ x{T)dT < 00 and x{t) G Loo if for 

each element of x{t), supj>g|xi(f)| < 00 , i = 1,..., m. 

2 Problem Statement 

Here we consider n physical agents described by double¬ 
integrator dynamics 

ii{t) =Vi{t), 

Vi{t)=Ui{t), i = (1) 

where Xi{t) G MP and Vi{t) G are, respectively, agent 
Es position and velocity, and Ui{t) is its control input. Let 

x{t) = [x{{t),...,xl{t)f andv{t) = 

An undirected graph G = (E, E) is used to character¬ 
ize the interaction topology among the agents, where 
V = {1,..., n} is the node set and E C V x V is the edge 
set. An edge {j, i) G E means that node i can obtain in¬ 
formation from node j and vice versa. Self edges {i, i) are 
not considered here. Let m denote the number of edges in 
E, where the edges (j, z) and (z,j) are counted only once. 
The adjacency matrix A = \aij\ G of the graph G is 

defined such that the edge weight = 1 if (j, i) G E and 
Oij = 0 otherwise. Eor an undirected graph, aij = aji. The 
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Laplacian matrix L = [lij] £ associated with A is 

defined as In = ^ij hj = where i ^ j. For 
an undirected graph, L is symmetric positive semi-definite. 
By arbitrarily assigning an orientation for the edges in G, let 
D = [dij] £ be the incidence matrix associated with 

G, where dij = — 1 if the edge ej leaves node i, dij = 1 
if it enters node i, and dij = 0 otherwise. The Laplacian 
matrix L is then given by L = DD^ [17]. 

Assumption 2.1 The undirected graph G is connected. 

Lemma 2.1 [17] Under Assumption 2.1, the Lapla¬ 

cian matrix L has a simple zero eigenvalue such that 
0 = Ai(L) < A 2 (L) < . •. < \n{L), where Ai(-) denotes 
the ith eigenvalue. Furthermore, for any vector y £ K" 
satisfying l^y = 0, X 2 {L)y’^y < y'^Ly < A„(L)y^y. 

Suppose that each agent has a time-varying input signal 
Tiit) £ i = 1,... ,n, satisfying 

r,{t) =<(f), 

ilit) =a:{t), ( 2 ) 

where (<) and a[ (t) £ are, respectively, the 

agent Fs input velocity and input acceleration. De¬ 
fine r{t) = [rl,...,rlf, v'^{t) = [vl^,... and 

..., Here we assume that the input 

signals are generated internally by the agents and each 
agent has access to its own input signal, input velocity, and 
input acceleration. 

We study the distributed average tracking problem for the 
double-integrator agents. The goal is to design Ui{t) for each 
agent to track the average of the input signals and input 
velocities, i.e.. 


1 X ^ 

^lim ||x*(f) - - '^rj{t)\\2 =0, 
oo n 

1 " 

lim \\v,{t) -= 0 , i = l,...,n, 

t—fOO Jl 

i=i 

where each agent has access to its own input information 
and has only local interaction with its neighbors. We are 
interested in controller design with reduced requirement on 
velocity measurements. First, we consider the scenario that 
each agent has communication capabilities but without the 
need for either absolute or relative velocity measurements. 
The agents achieves distributed average tracking for a wide 
class of input signals with bounded deviations among input 
accelerations in the presence of correct state initialization. 
Second, we consider the scenario that each agent has sensing 
but not necessarily communication capabilities without the 
need for relative velocity measurements. The agents achieve 
distributed average tracking even in the absence of correct 
state initialization. 


3 Velocity free distributed average tracking in the pres¬ 
ence of correct initialization 

In this section, we consider the case where the agents can 
initialize their initial conditions correctly. To remove both 
absolute and relative velocity measurements, we introduce 
the following filter for each agent 

Mt) = ~ 

j&Ni j&Ni 

jeNi 

Wi{t) = v^{t) - a ^ [xi{t) - Xj{t)], i = l,...,n, 

j&Ni 

(4) 

where Vi £ is an auxiliary filter variable, Wi £ is the 
filter output, and /3, a and 7 £ IR+ are control gains to be 
designed. 

We propose the following distributed control law for agent i 

u,(t) = - a - Wjit)] 

j&Ni j&Ni 

+ 7 X! sgnK(f) - tt;j(f)] 

j&Ni 

+ al(t), i = l,...,n. (5) 

For notational simplicity, we will remove the index t from 
variables in the reminder of the paper. 

Assumption 3.1 The deviations among input accelera¬ 
tions are bounded, i.e., ■sMPtg[o,oo)ll®r(0 ~®j(^)ll 2 < dj], 
i, j = 1,... ,n, where a]] £ M+. 

Assumption 3.2 


n n “ 

2=1 2 = 1 


Theorem 3.3 Using the control law given by (5), (4), and 
(3) for system (1), distributed average tracking is achieved 
asymptotically, provided that Assumptions 2.1, 3.1 and 3.2 
hold and the gains a, 7 and (3 are chosen such that 


ri 1 A2(L) -f 1 


^ A special choice is 2;i(0) = ri(0), Wi(0) = Vi(0) for all 
i = I,... ,n. 
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7 > (n - l)aS, 


1 + a^X^iL) 


^ 4A2(L)(aA2(L)-l)(a-l)’ 

where \ 2 {L) and Xn{L) are defined in Lemma 2.1. 


It can be proved that if p > and a > 1, then 


^l\ 2 {L) 1 , 1 r 1 

— - 1 1 > 0. Thus, using Shur com- 

1 ij “ [ij 

plement it is concluded that V is positive definite. The 
derivative of V along the trajectories of (7) and ( 8 ) is given 


Proof: The proof contains two steps. First, we prove that for 
each agent, Xi - ^ 'Ej^i ^ 0 and Wi - ^ Vj 0 . 

Then by showing that Xi — E”=i’"* ^ 0 and 

X]r=i ~ Sr=i concluded that 

- 7 Ej=i ^0 ^ Oandvi - ^ X;”=i fj ^ 0, i = 1,..., n, 
and hence distributed average tracking is achieved. 

Using the control law (5) for (1), we can get 


Xi =Vi, 

Vi=- a'^ [x^{t) - Xj{t)] + /? ^ [Wi{t) - Wj{t)] 
j&Ni j&N, 

+ 7 X! sgnKW i = l,...,n. 


(g) Ip] 

^np ^np 

V I I I £ 

* S ^np ^np ^np S 

Inp Inp (^Inp 

® ^p]^v + + £w^y ~ ® ^p\£x 

- a£y [L (g) Ip]£x + l^£w{^ ® Ip]£w 

+ iClD (g) /p]sgn([i7^ (g) Ip]£^) + + £^[L (g) Ip]£^ 

- a£'^[L (g) Ip]£y +C[L<Si Ip]£x - a^lL (g) Ip]£v 

- (^P£Z[L <S> Ip]£w - ajI.ZiD (g) /p]sgn([iA'^ (g) Ip]£yj) 

- a£'^a''- a‘^£'^[L(^ Ip]£y, (9) 

where [M (g) Ip]£x = [M"^ (g) /p]a; = [M (g) Ip]x = £x and 
using the same analysis [M^Ip]£y = £y and [M(g) = 

We can also analyze the term (1 — a)£'^a^ as 


Define £x = [M ^ Ip]x, £v = [M ^ Ip]v and = [M (g) Ip]w, 
where M = In - ^Inl^ and w = [wf,..., Wn]'^. Then 
the closed-loop dynamics ( 6 ) can be rewritten as 

£x —£vj 

£v=- a[L (g) Ip]£x + P[L (g) Ip]£w 

+ 'yID (g) /p]sgn([D^ (g) Ip]£w) + [M (g) 7p]a'', (7) 

where it can be easily proved that M xM = M, Lx M = L 
and D X M = D.lt follows from (3) and (4) that 

£nj=[L(g, Ip]£x - /3[T (g) Ip]£w 

- 7[£) (g) /p]sgn([i9^ (g) Ip]£w) - a'' - q;[T (g) /p]^„. 


Consider the following function 


(1 - =(1 - a)C[M^ (g) /p]c 


<{a - 1)\\[M (g) /p]Cu;|| 2 ||[AT (g) /p]a ’'||2 

i=l j=l,j^i 

( 1 A"U ^ ^ 

i=l j=lj^i 

n 

<(a-l)a;jmax{ V | |Cu.i - 112 } 

I 


i—1 jGNi 


i—1 jGNi 




fi[L 0 Ip] Inp Inp 

Inp Inp Inp 
Inp Inp ^Inp 


where /x € M+ and ^ . Since ^Jlnp = 0, 

using Lemma 2.1, we have 

1 r 

1 11 <S)iup)£ 

1 la 


where we have used the Assumption 3.1 and ||•|| 2 < ll'Hi 
to obtain second and last inequalities, respectively. 

Let fi = 2a — 1. Note that a > ^ 2 )Hl) ensures that p > 
Combining (9) and (10), we can get 

<£p£v + + (1 - a)iP[L ® Ip]£x 

+ {I - a)l3£l,[L ® Ip]£n, 

n 

- 7 (a - 1 ) 55 55 

i=l i&N, 
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i—1 j^Ni 

- aClL 0 Ip]^v - ® Ip]^v 

~^v ” 1 “ (1 ~ ^)ix [^ ® ^p\^x 

+ (l-a)/3^^[L(g)/p]Cu. 

- ^ “ C»i)’^sgn(^„i - 

*—1 j^Ni 

+ (“~if ~‘) »5 i:eii^--^-.ii. 

z —1 jGNi 

- aClL 0 /p]^« - 0 /p]C« 

Qll Qi2 Qi3 

Q 2 I <522 Q 23 C) ( 11 ) 

Q 3 I <532 Q 33 

where Qn = (1 - a)A 2 (L)/„p, Q 22 = [1 - Q;A 2 (L)]/„p, 

<523 = Q 32 = Inp — a^[i 0 /p], Q 33 = (1 ~ Oi)P\2{L)lnp 

Qi 2 = Qi 3 = Q 21 = Q 31 = 0 „p, and we have 
used 7 > (n — l)d^ and Lemma 2.1 for the last in¬ 
equality. We also used (g) Ip]sgn{[D'^ (g) Ip]^w) = 

Er=i T^jem to derive the condition of 

7 - 

Therefore, if the control gains a and (3 satisfy the constraint 
a > max{l, j^} and /3 > 4 A.(LKaA" ^-iKa-i) ^ it 
be seen that V is negative definite. Integrating both sides of 
(11), we can obtain € L 2 . Since V is positive 

definite and V is negative definite, ^v, G L^o- Thus, 

we can get ^xi ^vi ^ ^2 L L^o and ^x^ ^vi ^ ^ao- 
Using Barbalat’s lemma [18], it is concluded that ^x, 

—> 0 as f —)■ 00 which imply that Xi —^ 

Second, we analyze the terms X]r=i Vi- Con¬ 

sidering the closed-loop dynamics ( 6 ), under Assumption 
2 . 1 , the derivative of Sr=i 1 *® calcu¬ 

lated as 

n n 

n n 

i=l i=l 

where we have used the fact that the graph G is undi¬ 
rected. Using Assumption 3.2, it can be proved that 

i=i = T,i=i ri and Y^i^i = Ei=i <■ 

Combining the two parts shows that Xi —> ^ E "=1 Tj and 
7 Ej=i asymptotically. Therefore, distributed av¬ 
erage tracking is achieved asymptotically. ■ 


Remark 3.4 Compared with [13], using the algorithm de¬ 
fined by (3)-(5), the requirement for velocity measurements 
is removed. In addition, the only limitation on the input sig¬ 
nals is that the deviations among their accelerations are 
bounded. Therefore, the allowable input signals are more 
general than those in [13]. 

4 Distributed average tracking in the absence of neigh¬ 
bors’ velocity information and correct initialization 

In the previous section, the proposed algorithm solves the 
distributed average tracking problem in the presence of com¬ 
munication and accurate state initialization without velocity 
measurements. In this section we deal with the distributed 
average tracking problem without communication in the ab¬ 
sence of neighbors’ velocity information and correct initial¬ 
ization. 

We introduce the following filter 

Ti = ^ aij{xi - Xj) - fiwi - ysgn{wi) - nn 

j&N, 

- ( 12 ) 


Wi =Vi - a'^{Xi - Xj), i = l,...,n, (13) 

j&Ni 

where Vi G is an auxiliary filter variable, Wi G is the 
filter output, and k, a, y, /3 G are control gains. 

Assumption 4.1 Each input signal, its velocity and ac¬ 
celeration are bounded, i.e., sup^^j^ ^^\\ri(t )\\2 < f, 
*«^’te[ 0 ,oo)IK(i)ll 2 < and sup^^^Q^^•^\\dl{f ))\\2 < oT, 
i = 1,... ,n, where r, v’’, aT G K’*'. 

We propose the following control input for each agent as 

Ui=- n{xi - rf) - K{vi - vl) - a ^ (x* - Xj) j3w^ 

j&Ni 

-f 7 Sgn(wj)-f 4, i = l,...,n. (14) 

The following theorem presents sufficient conditions to solve 
distributed average tracking without neighbors’ velocity and 
accurate initialization. 

Theorem 4.2 Distributed average tracking is achieved for 
system (1) using (14), (13) and (12) asymptotically, under 
Assumptions 2.1 and 4.1 and provided that the gains k, a, 
7 and j3 are chosen such that 

K > 1, 


a > y/n, 
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7 > 


Q! + 1 , _ 

- + Kv +a j, 

a — yn 


^ 4 (a — 1 ) + (a — l)A2(i) 

(k - 1)^ + 2a^(K - 1)A„(L) + a'^A^(L) 

a\2{L) + K — 1 ’ 

where \ 2 {L) and \n{L) are defined in Lemma 2.1. 

Proof: Similar to Section III, the proof contains two 
parts. In first part, we prove that Xi — ^ 

~ n ^3 Wi —> 0 asymptotically. Define 

fx, ^v, M, and w as in the proof of Theorem 3.3. Now 
rewrite the dynamics (1) using (14) in vector form as 


iv = — Kfx + (g) Ip]r — + k[M (g) Ip]v'' 

- a[L (g) Ipl^x + P[M (g) Ip]w + 7 [M (g) Ip]sgn{w) 
+ [M(g)/p]a’'. 

Consider the following Lyapunov function candidate 



/il [L 0 Ip] -f /i 2 -^np Inp 

Inp Inp 

Inp Inp 




where /ii,p 2 € Using the same analysis as previous 
section, if 131 X 2 + P '2 > 1 and a > 1, U is positive definite. 
The derivative of V is given as 


V — pl^x ® Ipl^v + p2^x + '^ ^X 

~ ~ ^^x ® ^p]^x ~ k£,v ^X ~ 

— afy[L (g) Ipl^x - Kw'^^x + kw'^[M (g) Ip]r 
— Kw'^^y + kw'^[M (g) Ip]v^ + /3w'^[M (g) Ip]w 
+ 7 W^[M (g) /p]sgn('u;) + w'^[M (g) Ip]a^ 

+ ^2; [L (g) Ip]^x ~ Ol^x ® ^p\^v + [T <g) Ip^^x 

— [U <g) Ip]^v — ajIvFw — a'^vFs,gn{w) — anvFr 

— anyFy"^ — avFof — F'vF[L ® Ip\^y. (15) 


The term w '^{[M (g) Jp] (/tr + ku’’ + a’’ + 7 Sgn(w)) — a(^r + 
Ku’’ + o’" + 7Sgn(w))} can be analyzed as 


vF{[M (g) /p](Kr + KV^ + a’’ + 7 Sgn(rt;)) — a{Kr + kv'^ 
+ o’" + 7Sgn(u;))} 

<[(1 + a){Kf + KV^ + of) + 7-\/n] • ||w||2 — a7||t(;||i 

<[(1 + a)(ter + Ku’’ + a’’) + 7-\/n] • ||w||i — Q; 7 ||t(;||i. 

(16) 


Combining ( 15 ) and ( 16 ), if 7 > (ttr + kv’~ + a’’) and 
/ii = 2a — 1 and p2 = 2 k, we can get 

V <^y^v + 'uFfy - K^l^x + (1 - a)X 2 (LFx^x - 

— kvF^x — kvF + jdvFw — aX2{L)^y — ajdvFw 

- FnF[L (g) Ip]^y 

Pll P12 P13 ^x 

= F F P21 P22 P23 fv , ( 17 ) 

P31 P32 P33 W 

where Pn = [-k + (1 - a)A2(L)]/„p, P22 = [I - k - 
aA 2 (P)] 7 ^p, P12 — P2I — Onp P13 — -U31 — 2 ^ 7 np 7 

P23 = P32 = |[(1 — tv)/„p — FL], and P33 = ( 3{1 — a)Inp 
and we have used the fact that M — In is negative semi- 
definite and Lemma 2.1 for the inequality. 

If the control gains k, a, and /3 satisfy the constraints men¬ 
tioned in Theorem 4 . 2 , the matrix P and thus V are neg¬ 
ative definite. Integrating both side of ( 17 ), it is concluded 
that ^X7 and w € L2. Since V is positive definite and V 
is negative definite, ^x, and w G Lqo. Since ^1, ^2 and 
w are bounded, it is concluded that ^x, ^y, w —>■ 0 asymp¬ 
totically. Therefore it is proved that Xi — X 0 ’ 

n Sj=i d and Wi —)• 0 as f —)• 00. 


In second part, we prove that Er=iP and 

Er=i Er=i F asymptotically. Define the variables 


/~i \—^7T V— \Tl j X — \TL X — \Tl 

Si = T,^=l Xi - T,^=l + and 5*2 = Ei=i + “ Ei=i' 


we can get from (2) and (14) that 


=^ 2 , 

n n 

82 =-kS 2 -kSi+ j3'^Wi+-i'^sgn{wi). (18) 

i=l 

We then use input-to-state stability to analyze the system 
(18) by treating the term /3E”=i + tELi sgn(wi) as 

the input and Si and S 2 as the states. If k > 1, the matrix 

is Hurwitz. The system (18) with zero input 

is exponentially stable and hence input-to-state stable. Since 
tUi —>■ 0 as f —)■ 00 for each agent, it follows that S'! —> 0 
and S 2 -G 0, which implies that Er=i Er=i P 

and Er=i Er=i respectively. Employing the re¬ 
sult of the first part, it is concluded that ^ Ej=i 

and Vi ^ X Ej=i F and the distributed average tracking is 
achieved. ■ 



Remark 4.3 Note that the algorithm in [13] relies on cor¬ 
rect initialization and relative velocity measurements. Using 
the algorithm defined by (12)-(14), there is no need for ac¬ 
curate initialization and absolute instead of relative velocity 
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Fig. 1. The network topology is used to characterize the interaction 
among ten agents. 



measurements are sufficient. The corresponding trade-off is 
that this algorithm relies on the assumption of bounded in¬ 
put signals, input velocity and input acceleration. 


Fig. 2. The positions of 10 agents and the average of the input 
signals, where the states are initialized correctly. The initial average 
of input signals is represented as a dot while the initial positions 
of the agents are represented as squares. 


Remark 4.4 To implement the algorithm (3)-(5), each agent 
needs its own and neighbors’ positions and filter outputs as 
well as its own input acceleration. Hence communication is 
needed between neighbors, while there is no need for either 
absolute or relative velocity measurements. This algorithm 
needs correct initialization but can deal with more general 
input signals. To implement the algorithm (12)-(14), each 
agent needs its own position, velocity, and relative position 
between itself and neighbors, and its own input signal, in¬ 
put velocity and input acceleration. In this case, communi¬ 
cation is not required and the algorithm can rely on only 
local sensing, where there is no need for relative velocity 
measurements among neighbors. The second one does not 
use correct initialization but requires the input signal, its ve¬ 
locity and acceleration all to be bounded. As a result, both 
algorithms have their values depending on the application 
scenarios. Reduced requirement on velocity measurements 
is achieved in both algorithms to reduce costs. 


5 Simulation 

In this section, numerical simulation results are given 
to illustrate the effectiveness of the theoretical results 
obtained in Sections 3 and 4. It is assumed that there 
are ten agents (n = 10), where the network topol¬ 
ogy is described by Fig. 1. In the first case, the in¬ 
put acceleration for agent i is given by [0.1(sin(5f) 4- 
mod(f, 2)) X 0.1(cos(5f) 4- mod(t, 2)) x and the 
initial position and velocity of the agents are chosen as 
a;(0) = [0.1 X b, —0.2 x and u(0) = [0.2 x b, 0.1 x 6]^, 
where b = [-4,-3,-2,-1, 0,1,2, 3,4, 5]^. We let 

ri(0) = Xi{0) and t;[(0) = t;i(0), z = 1,..., 10. We denote 
the jth component of Xi as x^ . Similar notations are used 
for Vi, Ti, and v^. The control parameters for all agent are 
chosen as a = 20, 7 = 5 and /3 = 400. We simulate the 
algorithm defined by (3)-(5). Fig. 2 shows the positions of 
the agents and the average of the input signals. Clearly, all 
agents have tracked the average of the input signals in the 
absence of velocity measurements. Fig. 3 shows the veloc¬ 
ities of the agents and the average of the input velocities. 
We see that the distributed average tracking is achieved for 
the agents’ velocities too. 




Fig. 3. The velocities of 10 agents and the average of input signals, 
where the agents’ velocity are initialized correctly. 

In the second case, we simulate distributed average 
tracking when the agents’ positions and velocities are 
not initialized correctly. The input acceleration for each 
agent is described by [0.1 sin(f) x z, 0.1 cos(f) x z]^. 
We simulate the algorithm defined by (12)-(14). The 
initial values are set as x(0) = [0.1 x 6,—0.2 x 6]^ 
and z;(0) = [0.2 X 6, 0.1 X ri(0) = [0,0]^ and 

vl (0) = [—0.1 X z, —0.1 X iff, z = 1,..., 10, and the con¬ 
trol parameters are chosen for all agents as zt = 2, a = 10, 
7 = 50 and /3 = 450. Fig. 4 and Fig. 5 show that the dis¬ 
tributed average tracking is achieved for both agents’ po¬ 
sitions and velocities in the absence of neighbors’ velocity 
information and accurate initialization. 


6 CONCLUSIONS 

In this paper, distributed average tracking for a group of 
double-integrator agents has been studied. First a distributed 
discontinuous controller combined with a distributed bi¬ 
ter was proposed in the absence of velocity information to 
deal with input signals with bounded acceleration devia¬ 
tions, where the agents’ position and velocities must be ini¬ 
tialized correctly. Here only the relative position and biter 
output information were used in the control design besides 
each agent’s own position and input acceleration. The al¬ 
gorithm was then modibed to remove the requirement of 
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Fig. 4. The positions of 10 agents and the average of input signals 
in the absence of correct initialization. The initial average of input 
signals is represented as a dot while the initial positions of the 
agents are represented as squares. 




Fig. 5. The velocity of 10 agents and the average of input signals 
in the absence of correct initialization. 

communication and correct initialization under the assump¬ 
tion that the inputs, their velocities and their accelerations 
are bounded. In both cases, reduced requirement on velocity 
measurements was achieved and hence the cost for velocity 
measurements was reduced. 
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